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Phase change materials can be reversibly switched between amorphous and crystalline states and 
often show strong contrast in the optical and electrical properties of these two phases. They are 
now in widespread use for optical data storage, and their fast switching and a pronounced change 
of resistivity upon crystallization are also very attractive for nonvolatile electronic data storage. 
Nevertheless there are still several open questions regarding the electronic states and charge trans- 
port in these compounds. In this work we study electrical transport in thin metallic films of the 
disordered, crystalline phase change material GeiSb2Te4. We observe weak antilocalization and dis- 
order enhanced Coulomb interaction effects at low temperatures, and separate the contributions of 
these two phenomena to the temperature dependence of the resistivity, Hall effect, and magnetore- 
sistance. Strong spin-orbit scattering causes positive magnetoresistance at all temperatures, and a 
careful analysis of the low-field magnetoresistance allows us to extract the temperature dependent 
electron dephasing rate and study other scattering phenomena. We find electron dephasing due to 
inelastic electron-phonon scattering at higher temperatures, electron-electron scattering dephasing 
at intermediate temperatures, and a crossover to weak temperature dependence below 1 K. 



I. INTRODUCTION 

Phase change materials such as Ge 2 Sb 2 Te 5 or 
GeiSb2Te4 can be rapidly and reversibly switched be- 
tween the amorphous and crystalline state. This phase 
transformation is accompanied by a significant change of 
optical properties, which is exploited in rewritable optical 
data storage.^ Phase change materials are also attractive 
for non- volatile electronic memories, where a pron oun ced 
change of resistance 2 and fast electrical switching^! are 
advantageous. At present there are two challenges for 
commercial application of phase change materials in elec- 
tronic memories, and both are related to their electrical 
transport properties. First, the change of resistivity of 
the amorphous phase with time ("drift") is a disadvan- 
tage if multi-level storage concepts are to be realized. 
And second, the resistivity of the crystalline state in 
phase change materials is often low, so that high currents 
must be applied to heat the crystalline film to the melting 
point and subsequently bring the material to the amor- 
phous state. To identify phase change materials with a 
higher resistivity in the crystalline state, an in-depth un- 
derstanding of charge transport is a prerequisite. Recent 
transport studies^ have revealed a disorder-tuned metal- 
insulator transition in films of crystalline GeSbTe (GST) 
compounds. In this work it was shown that annealing 
of crystalline films to progressively higher temperatures 
is accompanied by a change from a negative temperature 
coefficient of resistivity (TCR) at low annealing tempera- 
ture, to a positive TCR at higher annealing temperature. 



This transition from an insulating to a metallic state is 
caused by a distinct increase of order upon annealing, 
and such an increase should have a pronounced impact on 
the low temperature resistance. In addition, a prominent 
change of the magnetoresistance at very low temperature 
is expected upon tuning the order, i.e. upon annealing 
of the samples. To investigate these predictions we have 
studied the low temperature electrical transport proper- 
ties of thin, quasi 2-dimensional (2D) films. 

Studies of electrical transport in low-dimensional dis- 
ordered conductors have revealed several novel phenom- 
ena that can appear at low temperature, including weak 
localization and antilocalization quantum interference 
(QI) and many body disorder enhanced electron-electron 
Coulomb interaction (EEI) effects^. The contributions of 
these two effects to the temperature dependence of the re- 
sistivity and Hall effect, as well as the magnetoresistance 
(MR), can be used to determine scattering mechanisms 
and other materials properties.^ In 2D films these con- 
tributions are logarithmic in temperature and have non- 
trivial dependences on the magnetic field. Magnetoresis- 
tance measurements in both perpendicular and parallel 
field geometries are useful for resolving these contribu- 
tions, since orbital QI effects are sensitive to field ori- 
entation and can be suppressed in an applied magnetic 
field, while the EEI effects are isotropic and generally 
much less sensitive to the magnetic field. These tech- 
niques are well established-^ but continue to be useful in 
the study of electrical transport in a variety of novel sys- 
tems, including grapheneP, unusual oxide heterointerface 
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structure d 9 * 10 ^, and topological insulator:^. 



There have been few experimental studies of low- 
temperature transport in phase change materials, and 
none investigating disorder-induced quantum effects in 
quasi-2D films, despite considerable theoretical and prac- 
tical interest in such results. Along with the recent stud- 
ies of metal-insulator transition in 100 nm thick films 5 , 
there have been limited investigations of longitudinal and 
Hall resistance in GST materials^ and also of weak lo- 
calization in the related compound GeMnTe^. 



In this article, we report the first study of disorder- 
induced quantum corrections, including weak antilocal- 
ization and enhanced EEI, in thin quasi-2D GST films. 
At high temperatures we observe a classical parabolic 
magnetoresistance and a large temperature independent 
contribution to the Drude resistivity. This contribution 
comes from static disorder and decreases monotonically 
with increasing annealing temperature. This large disor- 
der gives rise to a short electronic mean free path indica- 
tive of diffusive electronic conduction. At low tempera- 
tures (< 20 K), we find a resistance minimum followed 
by a small upturn in the resistance that is proportional 
to ln(T) as the temperature continues to decrease. The 
low temperature MR is positive everywhere and shows 
a sharp cusp that develops below 20 K, suggesting weak 
antilocalization QI arising from diffusive carrier trans- 
port in the presence of strong spin-orbit scattering. At 
high fields, we are able to recover the EEI contribution 
to both the resistivity and Hall effect. Finally, we an- 
alyze the MR measurements using established localiza- 
tion theories and determine contributions to electronic 
scattering. In future studies we also intend to investi- 
gate systematic changes of QI and EEI effects upon the 
transition from the metallic (weakly localized) to the in- 
sulating (strongly localized) state. The data and conclu- 
sions presented in this manuscript should help to under- 
stand the relationship between structural and electronic 
propertied^, recently reported high temperature magne- 
toresistive effects^, and the role of atomic vacancies in 
electronic propertied^. 



Our paper is organized as follows. In Sec. [IT] we review 
the experimental parameters for our GST films, and in 
Sec. |III| discuss the film dimensionality. Section |IV| con- 
tains an analysis of the low temperature resistance and 
Hall effect anomaly resulting from QI and EEI effects. 
Section [V] presents the measurements and analysis of the 
low temperature magnetoresistance in perpendicular and 
parallel field, and Sec. |Vl] contains a detailed discussion 
of the dominant scattering phenomena determined from 
the MR measurements. Finally, Sec. |VIl| summarizes our 
results and suggests open questions that may be of inter- 
est for future study. 
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FIG. 1: Upper panel: sheet resistance R xx versus temperature 
for three GeSbTe samples showing metallic slope, large residual 
resistance R(0), and a low-temperature minimum. Also shown are 
fits (solid black lines) to a Bloch-Griinciscn form (see text). Lower 
panel: Hall coefficient Rjj versus temperature for samples 1 and 2. 



II. EXPERIMENT AND SAMPLE 
CHARACTERIZATION 

Films of GeiSb2Te4 with thicknesses of 7.5 and 14 nm 
were deposited by sputtering onto Si substrates using sto- 
ichiometric GST targets, and capped with a 7 nm ZnS- 
Si02 layer (with a composition ratio of ZnS:Si02 equal 
to 80%:20%). The film thicknesses were controlled by 
adjusting the sputter time as confirmed using X-ray re- 
flectivity techniques. Films were annealed for 30 minutes 
in pure Argon gas flow at temperatures of 275 C (samples 
1 and 3) and 300 C (sample 2); 100 nm thick samples pre- 
pared identically were found to be in the hexagonal crys- 
tal phaseP After annealing, Hall-bar devices were pat- 
terned using conventional photolithography and Ar-ion 
milling techniques, with Ti/Au pads for electrical con- 
tact. The active area of the devices was 200 x 100 /im 2 . 

The longitudinal and Hall resistances were measured 
in several Quantum Design PPMS cryostats in applied 
magnetic fields of up to 9 T using standard four-point dc 
and low frequency lock-in techniques. Three completely 
independent experimental set-ups were used and all gave 
identical results. Care was taken to ensure that resis- 
tances were measured in the ohmic regime, especially at 
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TABLE I: Measured and calculated GST sample parameters. The sheet resistance R xx , carrier density n, mean free path £ e , 
and diffusion coefficient Dtr are measured or calculated from 15 K data. (We assume the sample 3 carrier density to be equal 
to sample 1.) The mobility extracted from the magnetoresistance (MR) fiMR and from the Hall effect HHaii were determined 
using 300 K MR data. Tmin is the resistance minimum temperature, and d the film thickness. 



Sample 


Tanneal 


d 




n 


4 


Tmin 




fJ-MR 


Dtr 




(°C) 


(nm) 




(10 20 cm- 3 ) 


(nm) 


(K) 


(cm 2 /V s) 


(cm 2 /V s) 


(cm 2 /s) 


1 


275 


7.5 


0.86 


1.9 


3.7 


14 


35.7 


61 


4.05 


2 


300 


7.5 


0.68 


2.2 


4.4 


12 


37.5 


73 


4.94 


3 


275 


14.0 


0.39 


(1.9) 


4.4 


9 


58.7 


62 


4.73 



low temperatures (< IK). Over 40 devices from the three 
films were fabricated and characterized at room temper- 
ature, and several for each film were further studied at 
low temperatures. All devices from the same film showed 
qualitatively identical behavior and only minor variation 
in their resistivity, Hall coefficient, and other parame- 
ters. All experimental data and formulas below report 
resistances R^j and conductances <Xy in 2D (sheet) val- 
ues. 

Figure [T] shows the measured sheet resistance versus 
temperature for all samples, which show a positive tem- 
perature coefficient of resistance and weakly metallic be- 
havior with a large zero-temperature residual resistance 
R(0). The decrease in the resistance between 300 and 
20 K is consistent with a reduction in electron-phonon 

scattering- Fig. [I] shows fits to the Bloch-Griineisen 

expression^ for the phonon contribution to electronic 
scattering, R(T)/R min — 1 ~ (T/9d) 5 at the lowest tem- 
perature. We find a Debye temperature 9d ~ 140 — 150 
K, consistent with independent measurements on compa- 
rable films^, and extract a value for the electron phonon 
coupling constant X ep ~ 0.1. 

The measured Hall resistance R xy is linear in the ap- 
plied magnetic field and is consistent with p-type charge 
carriers; we extract the Hall coefficient from high-field (B 
> 2 T) measurements of R xy . The lower portion of Fig. [11 
shows the Hall coefficient Rh = Rxy/B as a function of 
temperature for samples 1 and 2, which like the longitu- 
dinal resistance shows a low temperature minimum. The 
measured Hall coefficient at 15 K yields a carrier density 
n~ 2 x 10 20 cm~ 3 via \Rh\ = i^r t an d a Hall angle Tan 
Oh = Rxy/Rxx ~ 0.005 at a field of 1 T. 

The longitudinal conductance in the presence of a mag- 
netic field <r xx (B) can be calculated from R xx and R xy 
through 



cr X x{B) = 



R 2 



R 2 



R r 



(l + tan(6» H ) : 



(1) 



Since tan(##) -C 1, the conductance a xx pa 1/R XX , and 
the magnetoconductance Aa xx (B) — a xx (B) — a xx (0) pa 
^/Rxx{B) — 1/R XX (0). Below we plot and analyze the 
negative magnetoconductance, —Aa xx (B), which has the 
same sign as the MR: -Aa xx (B) sa AR XX (B)/R XX (Q) 2 . 

Now let us consider the normal state magnetoresis- 
tance. Figure [2] shows the measured MR at 300 K and 
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FIG. 2: Top panels: magnetoresistance for samples 1 and 2 at 
300 and 100 K. Lower panel: Kohler diagram showing scaling of all 
the magnetoresistance datasets in the upper panels. 



100 K for samples 1 and 2, defined as 

AR R XX (B) - R xx (0) 



R(0) 



Rxx(0) 



(2) 



At these temperatures, the MR is small, positive, and 
parabolic. Kohler's rule^ states that in a classical metal 
with one dominant scattering time, the magnetoresis- 
should be a universal function of the quan- 



tance 



AR 
R(0) 



tity w c r, where r is the transport scattering time, the 
cyclotron frequency uj c = eB/m* , and m* is the effective 
mass. At sufficiently low fields 



AR 

W) 



(3) 



consistent with the data presented in Fig. [2j From 
parabolic fits to the 300 K curve for sample 1 we de- 
rive a transport scattering time of r pa 1.4 x 10~ 14 s and 
a mobility fi pa 61 cm 2 /(V s), comparable to the cal- 
culated Hall mobility fiHall ~ 36 cm 2 /(V s) and some- 
what larger than that achieved in thick (100 nm) films 
annealed to the same temperatures^. For free electrons 
lu c t = BRft/Rxx and so the normal state MR should be 
a universal function of B/R xx ; such scaling for samples 
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1 and 2 is shown in the lower panel of Fig. [2] consis- 
tent with a free-electron, single scattering time picture 
of the classical electrical transport from 100 - 300 K. In 
our consideration of the low temperature MR below, this 
classical contribution was also apparent at the highest 
fields and was subtracted before proceeding with quanti- 
tative analysis. 

Based on the measured resistivity, Hall coefficient, and 
magnetoresistance, we calculated materials parameters 
including the mean free path i e and diffusion coefficient 
D tr as shown in Table [TT] for each sample, assuming a 
free-electron like picture with effective mass m* = 0.4 
m e and a valley degeneracy N v = 4P 



III. SAMPLE DIMENSIONALITY 

Before considering the WL analysis, let us discuss the 
dimensionality of our GST films. A thin film can be 
treated as two dimensional if the thickness d is smaller 
than the appropriate physical length scales. However, 
due to the many interactions and phenomena that are 
relevant in this analysis, there is no single measure of our 
film dimensionality. The electronic mean free path £ e for 
the samples is ~ 4 nm, less than d = 7.5 - 14 nm, and 
so the classical diffusive transport is three-dimensional. 
For QI effects, the relevant length scale is the dephas- 
ing length L^, related to the phase breaking time 
through Lcf, — U Dt^. In these samples the dephasing is 
dominated by inelastic scattering, both electron-phonon 
at high temperatures, and electron-electron scattering at 
lower temperatures. The rates for both of these scatter- 
ing process increase with temperature, and thus the 
decreases with temperature. At a sufficiently low temper- 
ature there should be a crossover to 2D behavior when 
L<f, ~ d, which (based on the analysis below) should occur 
at roughly 50-100 K. For EEI effects, th e relevant length 
is the thermal diffusion length, Lt = y Dh/kgT, which 
is equal to the film thicknesses at T ~ 40 K (sample 3) up 
to 140 K (sample 2). Thus despite being 3D with respect 
to classical transport, our films can be treated as 2D for 
QI and EEI phenomena below ~ 40 K. We restrict our 
analysis to this quasi-2D limit. 

The characteristic phonon wavelength \ p h 
hv s /kBT ~ 1 nm at 15 K, smaller than the film thick- 
nesses, indicating that the phonons cannot be considered 
two-dimensional and therefore constraining the theoret- 
ical predictions for the electron-phonon scattering rates. 
Finally, the characteristic magnetic length £b = \fh]AeB 
is equal to the 7.5 nm thickness of samples 1 and 2 at ~ 
3 T and equal to the sample 3 film thickness of 14 nm 
at ~ 1 T; at higher fields the diffusive transport is no 
longer in the quasi-2D limit. We now turn to an analysis 
of the anomalous upturn in the resistance apparent at 
temperatures below 20 K. 
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FIG. 3: Top panel: relative change in conductance versus tem- 
perature for all three GST samples along with linear fits in ln(T); 
the curves have been displaced vertically for clarity. Lower panel: 
conductance of sample 1 in zero field and applied magnetic fields of 
10, 20, 50, 100, and 200 mT. Inset: Slope of ln(T) fits for sample 1 
as a function of applied magnetic field on double logarithmic scales, 
showing saturation at fields above ~ 0.5 T. 



IV. WEAK LOCALIZATION AND LN(T) 
RESISTANCE 



A. Theoretical Background 

The low temperature electrical conductance a xx of a 
disordered metallic film in the presence of weak localiza- 
tion QI and EEI corrections can be written as 



a xx {B,T) = a + Act 



qi.eei 



(B,T) 



(4) 



= where cr = ne m l d is the (classical) Drude conductance, 
and the remaining term arises from disorder-induced QI 
and enhanced EEI effects. These corrections have been 
studied extensively^ and are especially pronounced in 
systems of reduced dimensionality. In 2D, both effects 
lead to a hi(T) contribution to the conductance. 

At high temperatures, when the electronic phase coher- 
ence time (t^) is very short, the diffusive charge trans- 
port in a disordered film can be treated classically. As 
the temperature is decreased, electron-phonon and other 
inelastic scattering events (which randomize the elec- 
tronic phase and thus limit r^) are suppressed. Once 
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70 becomes long relative to the elastic scattering time 
r e , quantum interference effects appear. One such ef- 
fect, weak localization (WL), results in an upturn in the 
resistance at low temperature arising from constructive 
interference of self-intersecting time-reversal-symmctric 
trajectories. Weak antilocalization (WAL) occurs when 
spin-orbit scattering is strong, and results in a decrease in 
the resistance. The WL correction, Aer(T) ~ In [r^/r e ], 
is sensitive to the dephasing time since only trajectories 
that are coherent within this timescale can contribute to 
the interference effect. The characteristic amplitude of 

2 

these effects is Go = ^rg , and if the inelastic scattering 
rate has a power law dependence on the temperature, 
t^ 1 ~ T p , then the WL correction is proportional to 
ln(T)P 

Aa QI (T) = GoptN vV In [T/T ] . (5) 

Here To is a reference temperature, N v is the number 
of degenerate valleys, rj is related to the strength of in- 
tervalley scattering, and the pref actor N v i] is of order 
unity. (N v r) — > 1 when intcrvallcy scattering is strong, 
and N v rj — > N v when intervalley scattering is weak.) 
Finally, t relates to the strength of spin-orbit scatter- 
ing, t = 1 in the weak spin-orbit scattering (WL) limit 
( T so 1 ^ T 4> )■> an d t = — | in the opposite (WAL) limit 

(^o 1 » V 1 )- 

The EEI correction to the conductance^^ also shows 
a ln(T) dependence 

Aa EEI (T) =G (l- ^F^j In [T/T ] (6) 

where F a is a 2D effective screening parame ter an d in the 
limit of strong spin-orbit scattering F a — > 0! 19 l 27 l Assum- 
ing that the QI and EEI effects are additive, the total 
zero-field contribution to the conductance is 

Aa(T) = G A\n[T/T a ] (7) 

where the prefactor A is given by 

A = N v r)tp + 1 - ^ (8) 

Modest applied magnetic fields suppress the WL correc- 
tion but not the EEI term, so we now consider our mea- 
surements of the conductance in zero and applied fields 
as a function of temperature. 

B. Temperature dependence of the conductance 

Figure [3] shows the conductance versus temperature for 
each sample. Also shown are linear fits in the tempera- 
ture range 0.7 K to 5 K, where the data are approximately 
linear in ln(T), with slopes A ~ 0.7. As shown in Eq. |8j 
these slopes contain contributions from QI and EEI ef- 
fects. Measuring this slope in applied magnetic fields 



allows separation of the two, since the QI effect can be 
suppressed with a moderate field, while the EEI contri- 
bution contributes only weakly to the MR and cannot be 
suppressed. The lower portion of Fig. [3] shows the con- 
ductance in applied fields ranging from up to 200 mT 
for sample 1, along with linear fits to the low-temperature 
data (dashed lines). The resulting logarithmic slopes for 
these fits (along with additional fits at higher fields, not 
shown on the plot) are shown in the inset. 

The slope of the longitudinal conductance is Al ow ~ 
0.70 in zero field, and saturates to Am ~ 1.2 at fields of 
roughly 0.5 T or more. Based on our analysis of mag- 
netoresistance measurements (see discussion below) we 
know that the characteristic dephasing field for sup- 
pression of the QI correction is less than 0.1 T throughout 
this temperature range; thus 0.5 T should be sufficient 
to suppress the QI localization effect. With the QI con- 
ductance correction slope A® 1 = Al ow - Am ~ -0.5, the 
remainder should be due to EEI: A EEI = A m ~ L2. 
The value for A® 1 is is consistent with strong interval- 
ley scattering (N v i] — > 1), t = (the strong spin-orbit 
scattering limit), and p w 1 as expected for electron- 
electron scattering dominated dephasing (see discussion 
below). In the limit of very strong screening and strong 
spin-orbit scattering F a — y 0, the EEI contribution A EEI 
should be 1. Our value of 1.2 indicates either that the 
parameter F a is negative or that there is an additional 
ln(T) contribution to the conductance appearing at low 
temperature beyond those considered here. 

Table |n] shows the collected results from the analysis of 
the zero-field and in-field resistance versus temperature 
measurements. We comment that the observed EEI ef- 
fects seen in these metallic samples showing WL behavior 
does not contradict the observations of Siegrist et alP, 
who highlighted the importance of disorder relative to 
electron interaction effects to explain the MIT in thicker 
films of this material. Since the QI and EEI corrections 
are expected to contribute to Rr in distinct ways, fur- 
ther consideration of temperature dependent Hall effect 
data can supplement the above analysis. 



TABLE II: Weak localization analysis of resistance versus 
temperature data for the three samples. The coefficient A 
is determined from zero-field data, and the EEI prefactor 
1 — jF a from data measured in appiied fields larger than ~ 
0.5 T. The QI prefactor tpN v rj is the difference between these 
values. 



Parameter 


Sample 1 


Sample 2 


Sample 3 


A 


0.71 


0.69 


0.56 


A eei = ! _ | F<r 


1.22 


1.18 


1.15 


A QI = t p N v rj 


-0.51 


-0.49 


-0.59 




FIG. 4: Relative change in the Hall coefficient Rh and longitu- 
dinal resistance at zero field (R X x(0)) and 5 T (R X x(5 T)) versus 
temperature for sample 1. Also shown are linear fits (dashed lines) 
to the low temperature region of each proportional to ln(T). Values 
for the slope S of the three linear fits are listed in the inset. 



C. Temperature dependence of the Hall effect 

Figure [4] shows the Hall coefficient Rh , along with the 
longitudinal resistance measured in zero field, R xx (0), 
and 5 T, R xx (5 T), all as a function of temperature. The 
three curves have been scaled by the resistance minimum 
value Rrnin'; plotted is: 



AR t _ Ri(T) — Ri 



R; 



Ri. 



(9) 



where Ri is either R xx (0), R XX (5T), or Rh . The behavior 
of Rh is similar to that analyzed for R xx in the previous 
section; Rh shows a minimum between 10 and 20 K fol- 
lowed by an upturn that is linear in ln(T). Also shown in 
Fig.|4]are linear fits to the low temperature ln(T) portion 
of each dataset; the slopes Si are listed in the inset. 

Both the QI and EEI localization effectiP^ contribute 
to the diagonal conductance a xx as described in Eq. |4j 
while only the QI effect contributes to the Hall conduc- 
tance a xv 



Aa QI 



(10) 



where here a E y is the classical Drude contribution to the 
Hall effect. When inverting the conductivity tensor cry, if 
only QI effects are present, then (to lowest order) there is 
no correction to Rh, while if only EEI effects are present, 
then the correction to the Hall coefficient Rh is twice that 
to the longitudinal resistance 



AR 



H 



Rh 



EEI 



f AR X 



(ii) 



Finally, if both QI and EEI corrections are present, then 
the Hall coefficient correction will include both QI and 



FIG. 5: Magnetoresistance for sample 1 (left) and 2 (right) from 
300 to 0.4 K; temperatures are indicated for each curve. 



EEI termfP 



AR 



Rh 



2A<j eei 



f AR X 



AaQ 1 + Aa EEI \ R x 



(12) 



Since our measurements of the Hall coefficient were 
performed at magnetic fields much larger than the de- 
phasing field (B ^> B^), we would expect this to suppress 
the QI term in Eq. [12] above, so that the correction to 
the Hall coefficient (ARh / Rh) should be 2 times the R xx 
correction as in Eq. [IT] The measured ratio (ARh / Rh) 
I (AR xx (hT)/R xx (VT)) = S H I S 5T is 1.4, close to the 
predicted value of 2 and confirming the presence of EEI 
corrections at high field. Studies of 2D electron systems 
in SP^and InZnO film:P!l observed ratios in the range 1-2, 
recovering the theoretical prediction of 2 in the R xx —> 
limit. 

In the low-field limit we expect from Eq. [12] to ob- 
serve (ARh/Rh) I (AR xx (0)/R xx (0)) = 2 S 5T /S = 3.6. 
The measured ratio Sh/Sq — 2.5 is comparable to this 
value, given that the Hall coefficient measurement was 
performed at high field. While it is difficult to obtain pre- 
cise Hall coefficient measurements in the zero- field limit, 
we did observe a qualitative increase in the slope Sh as 
B — >• 0. This increase would give a larger value of Sh/Sq, 
consistent with the theoretical prediction of Eq. [12] Ad- 
ditional study of the field and disorder dependence of the 
Hall coefficient at low temperatureshould further inform 
this analysis. Having resolved QI and EEI contributions 
to the resistance and Hall coefficient at low T, let us now 
consider the WL contributions to the magnetoresistance, 
with the final goal of understanding the scattering mech- 
anisms that govern the electronic dephasing processes in 
these materials. 
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V. MAGNETORESISTANCE ANALYSIS 

A. Perpendicular field 

Figure [5] shows the measured MR for samples 1 and 2 
at low temperatures, from 300 K to 0.4 K. In this temper- 
ature region, the Kohler rule scaling of the MR breaks 
down, and a pronounced cusp near zero field appears. 
Perpendicular magnetic fields suppress QI effects by de- 
stroying the phase coherence of time-reversal-symmetric 
paths. In the case of WAL this results in a positive 
MR, consistent with the data in Fig. [5] and indicating 
that we must consider spin-dependent scattering effects. 
The temperature and field dependent QI correction to 
the conductance, in the presence of spin-orbit and spin- 
flip scattering, was first described by Hikami, Larkin, and 
NagaokP (HLN) 



vqi{T,B) = -G a 



B J 2 \B J 2 \B 

(13) 

where ^(x) — tp(l/2 + x), ip(x) is the digamma function, 
a is a constant of order 1, and 61,2,3 are given by 




1E-3 



0.01 0.1 
^ CO 



FIG. 6: Negative magneto-conductance -Aa xx for sample 3 in 
perpendicular applied fields /jlqH ~ B^, at temperatures from 75 - 
1.8 K, plotted on a semi- logarithmic scale. The continuous curves 
show si ngle parameter fits to the theory of Hikami, Larkin, and 
NagaokipSI (HLN). The data and curves have been offset vertically 
for clarity. Inset: all magnetoconductance data from the main 
panel collapse to the theoretical HLN theory prediction (continuous 
curve) when plotted against ^lqH/ B^iT). 



B\ — B e + B so + B s 



Bi — ^B so 



-B s + Bi 



B-t — 2B S 



Bi 



(14) 
(15) 
(16) 



Here, B x = 4e j, r is the characteristic field correspond- 
ing to spin-orbit (so), spin-flip (s), elastic (e), and in- 
elastic (i) scattering processes. Note that the dephasing 
rate t^ 1 = 2t : T 1 + , and so we identify B3 above as 
the corresponding dephasing field, B^ — B 3 — 2B S + Bi. 
In the limit of strong spin-orbit and weak spin-flip scat- 
tering, B so » B Sl we may further simplify the above 
expression 



Aa(B) = -G a 



-A' 



1 B e 

2+5 

+ B d 



B 



B$ 
~B 



(17) 



Figure [6] shows the negative magnetoconductance de- 
termined for sample 3 at fields comparable to B^; note 
that the magnetic field is plotted on a logarithmic scale, 
and the curves are offset vertically. The magnitude of 
the magnetoconductance is of order of the WL prefactor 
Go- Figure [7] shows the negative magnetoconductance 
up to fields well above B^ for sample 1. Note that we 
have subtracted the (classical) parabolic contribution to 
the MR discussed in Sec. [Tl] above. At the highest fields 
measured, the sign of the MR remains positive though 
appearing to saturate. 

To reliably determine the characteristic fields appear- 



ing in Eq. 17 we first extract B^ from the low-field magne- 




FIG. 7: Negative magnetoconductance -Act versus applied per- 
pendicular magnetic field for sample 1 for temperatures from 50 K 
to 0.4 K, along with fits (continuous lines) to the WL theory of 
Hikami, Larkin, and Nagaoka as described in the text. 



high-field data. At low fields the classical magnetoresis- 
tance is negligible, and below B^ we expect ~ B 2 MR 
behavior due to QI, since in the limit of B <C {B so , B e ), 
Eq. [13] simplifies to 



y ' 48 I B A 



(18) 



toconductance, and then B so and a using the remaining 



Representative fits used to extract B^ using only the low- 
field portion (B < 5 x ZL,) of the magnetoconductance 
data are shown in Fig. [6] for sample 3. To determine 
the spin-orbit field we keep B^ constant from the low- 



field analysis and considered the full magnetic field range 
measured. The fits shown in Fig. [7] have only two free 
parameters - the spin-orbit field B so and the prefactor a 
appearing in Eq. [IT] (The elastic scattering field B e was 
fixed based on the estimated transport scattering time 
t « iCT 14 s, yielding B e ~ 10 T.) Deviations appear 
at fields of above ~ 3 T for this sample, however, this 
is already B e i ast i C /3 and therefore outside the region of 
validity for the perturbative (and diffusive-limit) HLN 
theory. In addition, above this field scale the film is no 
longer strictly 2D, since the magnetic length ig < d. The 
parameters B^ and B so from these fits are discussed be- 
low. High-field deviations from the HLN result are often 
attributed to spin-splitting or EEI contributions to the 
magnetoconductance. However, these effects should be 
isotropic and we observe no such phenomena appearing 
in the parallel field data considered next. 

B. Parallel fields 
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FIG. 8: Parallel field negative magnetoconductance -Aa xx 
for sample 1, along with fits (continuous lines) to the theory of 
Al'tshuler and Aronov. The curves have been vertically offset for 
clarity. 



When only considering orbital QI effects, a parallel 
magnetic field should not yield any MR with a strictly 
2D film. However, as first described by Al'tshuler and 
Aronov 21 (AA), a nonzero film thickness will allow for 
electron diffusion perpendicular to a field directed paral- 
lel to the film, and therefore lead to similar suppression 
of QI phenomena as in the perpendicular field case. Ac- 
cording to AA, the resulting magnetoconductance is 



Aa(B//) =G ln 



3l 

B-xBa 



(19) 



where B d = i|§ . When including spin-flip and spin-orbit 
scattering, this expression become^ 22 * 23 ! 



Aa(B f/ ) = G 



■In 



3i_ 

B2B d 



In 



3i_ 

B^B d 



(20) 

The magnetoconductance in parallel fields for sample 
1 was also measured for a similar set of temperatures and 
magnetic fields; these data are plotted in Fig. [8] If the 
measured parallel-field magnetoresistance arose from a 
small out-of-plane field appearing due to misalignment of 
the field orientation, then the data could be simply scaled 
to the perpendicular field results, but no such scaling is 
possible. While we cannot rule out the presence of a 
small out-of-plane contribution, the size of the measured 
magnetoconductance and lack of scaling indicates that 
this is not a significant contribution. Figure [8] also shows 
fits to the theory for QI contributions to the magneto- 
conductance in parallel field, Eq. 



20 with Ba and B s 



as free parameters. The resulting values for the dephas- 
ing field B^ are consistent with the perpendicular field 
analysis and are discussed below. In addition, the excel- 
lent agreement up to the highest magnetic fields suggests 
that isotropic effects, such as arising from EEI or Zeeman 
splitting, may be ignored in this field range. 
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FIG. 9: Characteristic dephasing and spin-orbit scattering 
B BO fields, extracted from fits to the magnetoconductance data, 
for all three samples in perpendicular (closed symbols) and parallel 
(open symbols) fields as a function of temperature. The upper 
panel shows B so on a linear scale, while the lower panel shows the 
strongly temperature dependent B^ on a logarithmic scale. 



Figure [9] shows the characteristic fields associated with 
dephasing and spin-orbit scattering, extracted from fits 
to the 2D WL theory as described above. The dephas- 
ing field B^ plotted on the lower panel shows reason- 
able agreement between the perpendicular and parallel 
field measurements for sample 1, and increases strongly 
with increasing temperature for all three samples. This 
strong temperature dependence is expected, since the 
field is proportional to the inelastic scattering rate which 
itself should increase as the temperature increases due to 
stronger electron-phonon and electron-electron scattering 
effects. We analyze this quantitatively in the next sec- 
tion. The spin-orbit scattering fields are roughly temper- 
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ature independent, particularly in the zero temperature 
limit where B so » B^, and appear to be consistent with 
established theory for surface scattering in the presence 
of high-Z atoms. 

Studies of interfacial spin-orbit scattering^ found that 
when surface scattering is the dominant elastic mecha- 
nism, B so can be expressed relative to the elastic scatter- 
ing field B e as B so — {af s Z) A B ei where Z is the atomic 
number and ctf 3 is the fine structure constant. Both 
antimony and tellurium have large Z; assuming stoichio- 
metric GeiSb2Te4 we estimate an average atomic num- 
ber Z ~ 49 and therefore with B e « 10 T would ex- 
pect B so ~ 0.16 T, somewhat smaller than the B so ~ 
0.5-2 T extracted above for the three samples. How- 
ever, the scaling of the spin-orbit times for samples 1 
and 2 (both 7.5 nm thick) yields T so ^ amp i e i /T so ^ amp i e2 ~ 
1.19, comparable to the ratio of elastic scattering times 
T e ,sampiei/Te, samp ie2 ~ 1.13 and consistent with the 
above modelsPSl 

The spin-orbit scattering rates are anisotropic and 
thickness-dependent, with Bjr Q ~ 0.8 T, and ML ~ 2.6 T 
for sample 1, and sample 3 (which is ~ 2 times thicker 
than 1 and 2) having Bjr D ~ 1.94 T. Similar studies of 
WAL in strong spin-orbit materials such as BP^ find no 
such anisotropy in B SOl and in contrast to our results that 
the spin-orbit scattering rate decreases with increasing 
thickness. Further study of high-field MR in both par- 
allel and perpendicular field may illuminate this discrep- 
ancy and inform ongoing study of spin-orbit coupling in 
phase change material compounds and in related topo- 
logical insulator materials such as Sb2Te 3 . 



VI. SCATTERING RATES 

Now let us consider the contributions to inelastic scat- 
tering in our samples, which we may compare with 
the characteristic dephasing rates r^" 1 extracted us- 
ing the magnetoresistance analysis above. In addition 
to temperature-independent dephasing due to spin-flip 
scattering from magnetic impurities or other extrin- 
sic effects®!, the two temperature-dependent contribu- 
tors to the dephasing are inelastic electron-phonon (c- 
p) and electron-electron (e-e) scattering. At sufficiently 
low temperatures, e-e scattering dephasing should dom- 
inate, since the e-p rate has a stronger temperature 
dependence^. 



A. Electron-electron scattering 

APtshuler, Aronov, and Khmelnitski (AAK)pl and 
other works^U found for the electron-electron scattering 
rate r,~ e _ e in the 2D limit (L T > d), the following ex- 
pression holds, 



,(e-e) 



(21) 



where R xx is the sheet resistance of the film, and the 
quantity G = 7 ^ zz . (We ignore an additional contri- 



butio 



to r. 



-l 

i,(e— e) 



that is proportional to T 2 , since it 



is smaller by a factor of ksT/Ep 



perature of 10 K, t 



i,(e— e) 



10 



11 



10~ 3 .) At a tem- 
, or comparable to 



the e-p dephasing rate (calculated below) at this tem- 
perature. At lower temperatures the dephasing field B^ 
scales with R xx for samples 1 and 2; this is expected since 
both the dephasing rate and the inverse diffusion coeffi- 
cient should be proportional to R xx . We therefore use the 
above theory to estimate the diffusion coefficient, since 



the only parameters in Eq. 21 are the film sheet resistance 
and fundamental constants. For sample 1, we obtain 7.2 
cm 2 /s, comparable to the 4.05 cm 2 /s estimated from the 
Drude scattering time. The temperature dependence of 
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the electron-electron dephasing rate is plotted in Fig 
as a broken line, and is the most significant contribution 
to the total rate at intermediate temperatures. This is 
consistent with the analysis of the resistance versus tem- 
perature data in Sec. Ill above, which showed dephasing 
linear in temperature between 2 and 10 K and a prefactor 
consistent with p = 1. Based on this consisten picture, 
we use the values for the diffusion coefficient extracted 
using the AAK theory Daak in the remainder of our 
analysis. 



B. Electron-phonon scattering 

There have been extensive measurements and analy- 
sis of e-p scattering contributions to dephasing in metal- 
lic and semiconducting thin films! 3 ^ In particular, vari- 
ous theories have proposed, and some experiments have 
observed, T 2 , T 3 , and T 4 power law dependence of the 
scattering rate. For example, Lawrence and Meador^l de- 
scribe the electron-phonon inelastic scattering rate T^~ e _ p 



in a dirty 2D film 



i,e— p 



T 



(22) 



where £(3) is the Riemann zeta function, and ujjj and 
0£> are the Debye frequency and temperature. With 0£> 
= 150 K and X ep = 0.1, at 10 K « 10 11 s" 1 , 

comparable to the dephasing rate in Fig. 10 Other the- 



ories predict power laws in temperature with exponents 
between 2-4 for electron-phonon scattering in a 2D film. 

In addition to the above T-linear and T x (with x ~ 2 
- 4) contributions to the electron dephasing rate, spin- 
flip scattering from magnetic impurities (or several other 
intrinsic or extrinsic processes, see discussion in e.gJ^J) 
may lead to a temperature-independent contribution to 
the dephasing. We therefore fit the measured dephasing 
rates, plotted in Fig. [lOj to the following form: 



'4> 



l (0) + C ee -T + C ep -T x 



(23) 



The zero-temperature rate 1 (0), and the T-linear coef- 
ficient C ee , were obtained from fits to the data below 5 K, 
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FIG. 10: Upper panel: dephasing scattering rate versus temper- 
ature on logarithmic scales for all three samples. The solid curves 
show fits to a theoretical model including e-p and e-e scattering and 
a phenomenological T-independent dephasing rate as described in 
the text. The data and fits for sample 2 (3) have been multiplied by 
a factor of 10 (100) for clarity. The dotted curves show the separate 
contributions for sample 1 from T-independent dephasing, T-linear 
electron-electron scattering, and T 2 6 electron-phonon scattering. 
Lower panel: low-temperature portion of the same data, plotted on 
linear axes. The continuous curves show the best-fit contributions 
from e-e and T-independent dephasing terms. 



the dephasing rate or a temperature dependence that is 
weaker than the AAK T-linear theory. Such saturation 
or weak T-dependence has been observed at low temper- 
atures in many studies^, and various mechanisms have 
been considered. Typically this effect is described by a 
temperature-independent contribution to the dephasing 
rate, r < ^" 1 (0), although the physical relevance of such a 
description in the limit of zero temperature may be un- 
clear. Saturation effects may include both intrinsic mech- 
anisms such as noise from two- level systems^ or spin- flip 
scattering, and extrinsic ones, such as an ambient mag- 
netic field, electron heating, or microwave noise effect™ 
Lin, Li, and Zhong observed a systematic dependence of 
T 1 (0) 011 D in disordered metals; based on this trend and 
with D w 5 - 10 cm 2 /s the predicted saturation rate is 
k. 2xlO n -2xl0 12 s~ 1 , considerably higher thanrT^O) 



3 - 4 x 10 



9, 



Recent studies of Sn-doped In 2 03_ 



films^S report a 1 (0) that decreases with increasing dis- 
order, whereas sample 3 above (which has the smallest 
value of R xx ) shows no sign sign of a temperature inde- 
pendent dephasing contribution. The magnetic field cor- 
responding to our tT (0) is ~ 10 Gauss, and we estimate 
that residual fields within our cryostats are much less 
than this value. If the zero-temperature dephasing rate 
is due to spin-flip scattering, then this implies a magnetic 
impurity concentration of roughly 10 ppirP^ (See also the 
discussion in McGinnis and ChaikerpZl). 

Finally, we comment on the departure from ln(T) be- 
havior in our K xx versus T data at the lowest tempera- 
tures, visible below 0.5 K in Fig. [3] for samples 1 and 2. 
A dephasing mechanism such as spin-flip scattering will 
suppress the QI correction but not affect the EEI contri- 
bution. In such a case the slope A would become larger 
since the QI and EEI ln(T) contributions are of opposite 
sign. However, we observe a downturn in the resistance 
(upturn in a xx ) in the T— > limit, suggesting that the 



while the parameters x and C ep for the electron-phonon 
contribution were obtained using only the high temper- 
ature (T > 10 K) data. The resulting coefficients, along 
with other analysis results, are collected in Table [ill] and 
the individual contributions for sample 1 are shown as 
dotted curves in Fig. [Tol 



C. Dephasing in the T — > limit 

The lower panel of Fig. [10] plots the low temperature 
dephasing rate data for all three samples along with the 
theoretical T-linear e-e dephasing and temperature in- 
dependent contributions. In an ideal metallic film e-e 
scattering should provide the only mechanism for elec- 
tronic dephasing at low temperature, and so the rate 
should continue to decrease as the length scale over which 
interference effects occur increases. Any cutoff in this 
length scale, such as that due to finite sample size or 
another scattering mechanism, will lead to saturation of 



TABLE III: Weak localization magnetoconductance analysis 
results. Dtr and Daak are the diffusion coefficient as deter- 
mined from a free electron picture and by scaling the mea- 
sured dephasing rate to the Al'tshuler-Aronov-Khmelnitski 
prediction (see text). B e and B so are the elastic and spin-orbit 
scattering transport fields, and t so is the spin-orbit time. The 
parameters 7"7 (0), C ee , C ep , and x describe the temperature 



dependent dephasing rate (see Eq. 23 in the text 



Parameter 



(units) 



Dtr 

Daak 
B e 

i 3 a 

C, 



X 



(cm 2 /s) 
(cm 2 /s) 
(T) 
(T) 
(T) 
(ltr 13 s) 
(10 9 s" 1 ) 
(10 10 s^/K) 
(10 7 s" 1 /^") 



Sample 1 Sample 2 Sample 3 



4.05 


4.94 


4.73 


7.0 


10.1 


11.1 


20.0 


12.6 


11.1 


0.81 


0.69 


1.94 


2.6 






2.9 


2.36 


0.76 


2.7 


4.8 





1.18 


1.01 


0.70 


5.9 


17.3 


10.5 


2.78 


2.56 


2.78 
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electronic system itself is unable to cool effectively at 
these temperatures. Assuming that the ln(T) behavior 
does remain, we scaled the above scattering rates with 
'effective' values for the temperature assuming continued 
ln(T) behavior. While the tendency towards saturation 
becomes less pronounced after this analysis, the temper- 
ature dependence is still weaker than T-linear, and we 
conclude that insufficient cooling is not responsible for 
the T-independent contribution to the dephasing t^~ 1 (0) 
in the above analysis. 

VII. SUMMARY AND CONCLUSIONS 

Our analysis has shown that even the most metallic 
phase change films studied here are governed by strong 
defect scattering and hence qualify as 'dirty metals', with 
a rich physics of QI and EEI effects. Significantly, in 
these materials we have a simple way of tuning the dis- 
order without changing composition. This provides an 
excellent opportunity to study systematic trends as a 
function of increasing disorder. By analyzing three sam- 
ples with different room temperature sheet resistances we 
are able to study disorder and thickness dependence and 
anisotropy of spin-orbit scattering, strengthening of the 
electron-phonon scattering rate with increasing disorder, 
and suggestions of a sub-T-linear power law governing 
electronic dephasing in the low-temperature limit. Con- 
tinuing such an analysis of more disordered films towards 
the limit of strong localization promises new insight into 
disorder induced MITs. In particular we expect the weak 
localization approximation to break down. Our prelimi- 
nary data analysis already reveals two surprising facets, 
a departure from ln(T) behavior and apparent saturation 
in the zero-temperature limit. 

We have studied the low temperature magnetotrans- 



port properties of disordered thin films of the phase- 
change compound GeiSb2Te4 that have been annealed 
to be weakly metallic. We observe clear signatures of 
weak antilocalization and disorder enhanced electron- 
electron interaction effects in the resistivity, magneto- 
conductance, and Hall effect measurements. Using estab- 
lished WL theory, we are able to extract several impor- 
tant materials parameters, including characteristic spin- 
orbit and inelastic scattering rates. We observe a ln(T) 
quantum correction in Rh 1.4 times larger than that seen 
in R XX1 comparable to the factor of 2 expected when 
the QI effects are suppressed at high field. The inelas- 
tic scattering is dominated by phonons at high temper- 
atures, electron-electron scattering at low temperature, 
and shows a weak temperature dependence at the lowest 
temperatures studied. The existence of saturation or a 
cutoff in the dephasing length at very low temperature 
has been studied extensively, and both intrinsic (spin- 
flip scattering due to magnetic impurities) and extrinsic 
(coupling to external dissipative phenomena) sources are 
often considered. Given the unresolved nature of dephas- 
ing mechanisms at low temperature, and the relevance of 
these issues to fundamental questions about 2D metal- 
lic phenomena, further study of these issues may prove 
fruitful. 
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